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Chapter  1 


INTRODUCTION 

The  swash  zone  on  a  beach  forms  the  boundary  zone  between  the  surf  and 
backshore.  A  predictive  model  for  swash  dynamics  is  required  for  modeling  the 
erosion  and  recovery  of  natural  and  artificial  beaches.  Swash  data  on  natural  beaches 
provided  the  qualitative  understanding  of  swash  dynamics  [e.g.,  Guza  and  Thornton 
(1982);  Holman  and  Sallenger  (1985);  Douglas  (1992)]  but  it  may  not  be  possible  to 
develop  simple  empirical  relationships  in  view  of  the  various  possible  combinations 
of  different  beach  profiles  and  incident  wave  characteristics.  On  the  other  hand, 
the  time- dependent  numerical  model  based  on  the  finite-amplitude  shallow-water 
equations  including  bottom  friction  (Kobayashi  and  Wurjanto  1992)  has  been  shown 
to  be  capable  of  predicting  the  swash  characteristics  on  a  natural  beach  within  errors 
of  about  20%  (Raubenheimer  et  al.  1995).  However,  this  shallow- water  model  can 
not  be  extended  to  intermediate  water  depth.  Furthermore,  the  time-dependent 
numerical  model  requires  significant  computation  time  to  resolve  the  breaking  wave 
profiles  varying  rapidly  in  time  and  space.  This  leads  to  the  study  of  possible 
statistical  methods  for  predicting  swash  dynamics. 

1.1  Previous  Statistical  Work 

The  time-averaged  models  for  random  waves  represented  by  the  root-mean- 
square  wave  height  (Battjes  and  Janssen  1978;  Thornton  and  Guza  1983)  or  ex¬ 
pressed  as  the  superposition  of  regular  waves  (Dally  and  Dean  1986;  Mase  and 
Kobayashi  1991)  are  more  efficient  computationally  at  the  expense  of  the  loss  of 


detailed  temporal  information,  such  as  the  skewness  of  the  wave  profile,  which  is 
regarded  to  be  an  important  characteristic  in  cross-shore  sediment  transport  [e.g., 
Guza  and  Thornton  (1985)].  Moreover,  the  existing  time-averaged  models  may  con¬ 
siderably  underpredict  the  wave  setup  and  root-mean-square  wave  height  in  the 
swash  zone  (Cox  et  al.  1994). 

Guza  and  Thornton  (1985)  computed  the  various  moments  of  the  fluid  ve¬ 
locity  field  measured  on  a  gently  sloping  beach  in  estimating  nearshore  sediment 
transport  rates.  The  time  averaging  and  probabilistic  averaging  of  a  random  vari¬ 
able  were  assumed  to  be  equivalent.  The  comparisons  of  the  computed  moments 
with  monochromatic  and  linear  (Gaussian)  random  wave  models  indicated  major 
shortcomings  in  the  monochromatic  representation  of  random  waves.  The  Gaus¬ 
sian  model,  with  given  mean  and  standard  deviation,  was  shown  to  predict  the 
even  moments  fairly  accurately  but  to  be  incapable  of  predicting  the  odd  moments 
associated  with  nonlinearities.  In  relatively  deep  water,  attempts  were  made  to 
predict  the  skewed  distribution  of  the  free  surface  elevation  using  higher-order  ex¬ 
pansions  [e.g.,  Longuet-Higgins  (1963);  Huang  and  Long  (1980);  Ta5Tun  (1980)].  In 
shallow  water,  Bitner  (1980)  compared  a  Gram-Chalier  expansion  with  field  data. 
This  higher-order  distribution,  with  measured  skewness  and  kurtosis,  predicted  the 
measured  free  surface  probability  density  functions  only  slightly  better  than  the 
Gaussian  distribution. 

1.2  Objective 

The  probabilistic  approach  has  recently  been  neglected,  probably  for  lack 
of  a  reasonably  simple  skewed  distribution,  in  favor  of  time-dependent  numerical 
modeling.  In  this  study,  the  exponential  gamma  distribution  [e.g..  Gran  (1992)]  with 
given  mean,  standard  deviation,  and  skewness  is  adopted  to  describe  the  probability 
density  function  of  the  free  surface  elevation  from  outside  the  surf  zone  through  the 
swash  zone.  The  range  of  the  skewness,  s,  allowed  in  this  distribution  is  s  =  0  —  2. 


The  skewness  determines  the  shape  of  this  distribution  and  the  , distribution  becomes 
Gaussian  for  s  =  0  and  exponential  for  s  =  2.  The  Gaussian  distribution  may 
be  appropriate  in  relatively  deep  water,  far  seaward  of  the  breaker  zone,  where 
local  nonlinearity  may  be  negligible.  The  exponential  distribution  may  be  a  good 
approximation  in  the  swash  zone  where  the  lower  limit  of  the  free  surface  elevation 
is  limited  by  the  bottom  elevation.  Since  there  is  no  theoretical  justification  for  the 
adopted  exponential  gamma  distribution,  an  experiment  was  conducted  to  evaluate 
the  applicability  of  the  exponential  gamma  distribution,  as  well  as  to  examine  the 
cross-shore  variations  of  the  mean,  standard  deviation,  and  skewness  of  the  free 
surface  elevation  from  outside  the  surf  zone  to  the  swash  zone. 

In  the  following,  the  exponential  gamma  distribution  is  explained  first  be¬ 
cause  the  present  experiment  was  designed  to  assess  the  cross-shore  variation  of  this 
distribution.  The  laboratory  experiment  is  then  described  concisely.  The  measured 
free  surface  elevation  and  velocities  are  analyzed  to  obtain  the  probability  distribu¬ 
tions  of  these  variables.  A  summary  of  this  study  has  been  submitted  for  publication 
(Kobayashi  et  al.  1996). 
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Chapter  2 


EXPONENTIAL  GAMMA  DISTRIBUTION 


®  The  free  surface  elevation,  rj,  above  the  still  water  level  (SWL)  is  assumed 

to  be  a  stochastic  process  satisfying  the  conditions  of  stationality  and  ergodicity 
[e.g.,  Goda  (1985)].  The  time  and  probabilistic  averaging  of  the  stochastic  vari- 
^  able,  denoted  by  an  overbar  in  the  following,  can  then  be  regarded  to  be  equivalent. 

The  probability  density  function  of  rj  is  assumed  to  be  expressed  by  the  exponen¬ 
tial  gamma  distribution  [e.g..  Gran  (1992)].  The  mean,  standard  deviation,  and 
skewness  coefficient  of  t]  are  denoted  hy  f],  a  and  s,  respectively.  It  is  assumed  that 
0  <  s  <  2  on  the  basis  of  available  field  and  laboratory  data  (Bitner  1980;  Thompson 
1980;  Huang  and  Long  1980;  Goda  1985;  Mase  and  Kobayashi  1991;  Raubenheimer 
et  al.  1995). 


2.1  Exponential  Gamma  Distribution 

The  normalized  free  surface  elevation,  j/,,  is  defined  as  r)^  =  {t}  —  fj)/a^  whose 
mean  and  standard  deviation  are  zero  and  unity  respectively.  The  skewness,  s,  of 
7/.  is  the  same  as  s  of  t].  The  exponential  gamma  distribution  for  t;,  can  be  shown 
to  be  expressed  as 

fin.)  =  lr(c.)l-'^/W  (-£-»)  (2.1) 

with 

y  =  7?.  -  tp{a)  (2.2) 
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in  which  a  =  shape  parameter;  F  =  gamma  function;  V’  =  digamma  function;  and 
Ip'  =  trigamma  function.  The  relationship  between  s  and  a  is  given  by 

s  =  -rw  [-/>'(«)i'‘'‘  (2-3) 

in  which  xp"  =  tetragamma  function.  The  kurtosis  K  of  77,  and  rj  is  given  by 

K  =  r{a)  +  3  (2.4) 

in  which  ip'"  =  pentagamma  function. 

The  gamma  and  related  functions  are  explained  in  Abramowitz  and  Stegun 
(1972)  and  tabulated  by  Gran  (1992).  Moreover,  the  values  of  r(a),  ip{a),  ip'{a), 
ip"{a),  and  ip'"{a)  for  given  a  can  be  found  using  Mathematica  (Wolfram  1991).  The 
value  of  a  for  given  s  can  be  obtained  by  solving  (2.3)  using  an  iteration  method 
(Gran  1992).  Figure  2.1  shows  the  variations  of  a  and  {K  -  3)  with  respect  to  s 
where  use  can  be  made  of  a  =  (0.5  +  s~^)  for  s  <  0.4  (Gran  1992).  The  shape 
parameter,  a,  decreases  from  infinity  at  s  =  0  to  zero  at  s  =  2.  The  kurtosis  K 
increases  from  /F  =  3  at  s  =  0  to  =  9  at  s  =  2. 

2.2  Limits  of  the  Exponential  Gamma  Distribution 

For  5  =  0,  equation  (2.1)  reduces  to  the  Gaussian  distribution 

On  the  other  hand,  for  s  =  2  equation  (2.1)  becomes  the  exponential  distribution 

for  5  =  2  (2.6) 

which  is  limited  for  the  range  rj:,  =  [{rj  -  rj)la]  >  -1.  This  lower  limit  may  be 
interpreted  in  terms  of  the  instantaneous  water  depth  h  =  {rj  —  z\P)  where  z\,  =  fixed 
bottom  elevation  which  is  taken  to  be  positive  above  SWL.  It  is  required  that  h  >0 
and  hence  rj  >  Z},.  The  lower  limit  of  77  for  s  =  2  may  be  assumed  to  be  imposed  by 


Figure  2.1:  Shape  Parameter,  a,  and  Kurtosis,  /i ,  as  a  Function  of  Skewness,  s 

Zb,  resulting  in  [fj  —  a)  =  Zb.  Since  the  mean  water  depth  h  is  given  hy  h  =  {fj  —  Zb), 
this  intuitive  assumption  requires 

h  =  a  for  s  =  2  (2-7) 

2.3  Limitations  of  Available  Data 

One  of  the  objectives  in  this  study  is  to  investigate  whether  equations  (2.6) 
and  (2.7)  hold  in  the  swash  zone.  The  available  surf  zone  laboratory  data  [e.g., 
Mase  and  Kobayashi  (1991)]  indicates  a  landward  decrease  in  skewness  s  toward 
the  still  water  shoreline.  It  is  hence  uncertain  whether  s  increases  landward  of  the 
still  water  shoreline,  and  approaches  two  in  the  swash  zone  as  implied  by  (2.6). 
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Moreover,  Raubenheimer  et  al.  (1995)  measured  the  free  surface  oscillations  using 
runup  wires  placed  parallel  to  a  natural  beach  face  and  obtained  small  and  negative 
values  of  s  for  the  sea  swell  frequency  range.  However,  the  free  surface  elevation  t] 
used  herein  corresponds  to  that  measured  by  a  vertical  gage  with  its  lower  limit  on 
the  beach  face. 

The  root-mean-square  wave  height,  is  defined  in  this  paper  using  the 

standard  deviation,  <t,  (Thornton  and  Guza  1983;  Battjes  and  Stive  1985) 

iJrms  =  \/8a  (2.8) 

If  equation  (2.7)  holds  in  the  swash  zone,  Hymslh  =  \/8  which  is  much  larger  than 
Hrmsih  0.42  in  the  inner  surf  zone  measured  on  a  natural  beach  by  Thornton  and 
Guza  (1982).  On  the  other  hand,  the  free  surface  elevation  time  series  measured 
by  vertical  capacitance  wave  gages  buried  partially  in  the  sand  beach  in  the  SU¬ 
PERTANK  Project  (Kriebel  1994)  indicated  that  Hmo/v  =  3  -  4  at  the  still  water 
shoreline,  where  Hmo  —  zero-moment  wave  height  given  by  Hmo  =  4o- .  Since  =  0 
and  h  =  7]  aX  the  still  water  shoreline,  this  observed  range  corresponds  to  i/rms/ h  = 
2.1  -  2.8.  This  suggests  that  equation  (2.7)  may  be  reasonable  in  the  swash  zone, 
and  that  the  value  of  Hrms/h  may  increase  appreciably  from  the  inner  surf  zone  to 
the  swash  zone. 


Chapter  3 


EXPERIMENTAL  PROCEDURES 

Experiments  were  conducted  to  obtain  a  dense  array  of  measurements  for 
free  surface  elevation  and  mid-depth  horizontal  fluid  velocity  in  the  shoaling,  surf, 
and  swash  zones. 

3.1  Experimental  Setup 

The  experiments  were  conducted  in  a  wave  tank  that  was  30  m  long,  2.44  m 
wide,  and  1.5  m  high  with  a  constant  water  depth  of  76.2  cm  (30  in).  Irregular 
waves,  based  on  the  TMA  spectrum  (Bouws  et  al.  1985)  using  linear  wave  theory 
and  random  phases,  were  generated  with  a  piston- type  wave  paddle.  A  rock  beach 
with  a  1:8  slope  was  located  opposite  the  wave  paddle  to  absorb  waves  propagating 
landward.  A  divider  wall  was  constructed  along  the  center  line  in  the  tank  to  reduce 
the  effects  of  re-reflection  from  the  wave  paddle.  The  experimental  setup  is  shown 
in  Figure  3.1.  A  plywood  beach  with  a  1:16  slope  was  installed  in  the  1.22-m-wide 
flume.  Eight  capacitance  wave  gages  were  used  to  measure  the  temporal  variations 
of  free  surface  elevations.  One  acoustic-Doppler  velocimeter  (Kraus  et  al.  1994)  was 
employed  to  measure  the  temporal  variations  of  fluid  velocities.  The  sampling  rate 
for  the  free  surface  and  fluid  velocity  measurements  was  20  Hz. 

Three  sets  of  tests  were  conducted  on  the  1:16  slope.  The  wave  paddle  motion, 
controlled  by  a  computer,  was  identical  for  each  test.  Wave  gages  1,  2,  and  3  were 
fixed  for  each  test  and  used  to  separate  the  incident  and  reflected  waves  using  the 
procedure  described  by  Kobayashi  and  Raichle  (1994).  Wave  gage  1  was  located  at 
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Figure  3.1:  Wave  Tank  Experimental  Setup 


the  water  depth  below  SWL,  d  =  75  cm,  for  tests  1  and  2  with  the  incident  wave 
spectral  peak  period  Tp  =  1.5  and  2.8s,  respectively.  The  horizontal  coordinate  x 
in  this  paper  is  taken  to  be  positive  landward  with  re  =  0  at  d  =  75  cm.  To  better 
resolve  the  incident  and  reflected  waves  for  test  3,  with  Tp  =  4.7  s,  wave  gage  1  was 
moved  to  rc  =  — 1  m  where  d  =  76.2  cm  above  the  horizontal  bottom.  For  all  three 
tests,  wave  gages  2  and  3  were  located  at  rc  =  0.38  m  with  d  =  72.6  cm  and  at  rc  = 
0.50  m  with  d  =  71.9  cm  respectively. 


3.2  Free  Surface  Tests 


Each  test  set  consisted  of  17  runs  performed  for  the  free  surface  measurements 
along  the  center  line  on  the  1:16  slope.  The  duration  of  each  run  was  400  s,  of  which 
the  initial  transient  duration  of  75  s  was  removed.  The  remaining  325  s  of  each  time 
series  was  used  for  the  subsequent  data  analyses.  The  runs  from  each  test  set  were 
numbered  1  through  17  for  the  free  surface  tests.  For  wave  gages  4-6,  the  horizontal 
spacing  of  the  two  adjacent  gages  was  0.16  m,  corresponding  to  a  1-cm  increment  of 
the  water  depth,  d.  These  three  gages  were  moved  in  the  region  of  x  =  1.44-9.44  m, 
where  d  =  16-66  cm,  to  measure  the  free  surface  oscillations  at  a  horizontal  interval 
of  0.16  m.  Simultaneously,  wave  gages  7  and  8  were  partially  immersed  in  gage 


Table  3.1:  Depth  of  Water  at  each  Gage  Position  for  Free  Surface  Test  1 


Run 

Water  Depth  d  (cm) 

Bottom  Elevation  Zi,(cm) 

Gage  4 

Gage  5 

Gage  6 

Gage  7 

Gage  8 

1-1 

66 

65 

64 

-12.5 

-10.0 

1-2 

63 

62 

61 

-12.5 

-10.0 

1-3 

60 

59 

58 

-10.0 

-7.5 

1-4 

57 

56 

55 

-10.0 

-7.5 

1-5 

54 

53 

52 

-7.5 

-5.0 

1-6 

51 

50 

49 

-7.5 

-5.0 

1-7 

48 

47 

46 

-5.0 

-2.5 

1-8 

45 

44 

43 

-5.0 

-2.5 

1-9 

42 

41 

40 

-2.5 

0.0 

1-10 

39 

38 

37 

-2.5 

0.0 

1-11 

36 

35 

34 

0.0 

2.5 

1-12 

33 

32 

31 

0.0  n 

2.5 

1-13 

30 

29 

28 

2.5 

5.0 

1-14 

27 

26 

25 

2.5 

5.0 

1-15 

24 

23 

22 

5.0 

7.5 

1-16 

21 

20 

19 

5.0 

7.5 

1-17 

18 

17 

16 

5.0 

7.5 

wells  as  illustrated  in  Figure  3.1.  The  horizontal  spacing  of  the  16  gage  wells  was 


0.4  m  in  the  region  10  m  <  a;  <  16  m,  corresponding  to  the  2.5-cm  increment  of  the 
water  depth  d.  Tables  3.1  -  3.3  show  the  water  depth  at  each  gage  position  for  tests 
1-3,  respectively.  Note  that  the  water  depth  is  given  by  d  for  gages  4-6  and  by  Zb, 
where  Zb  =  —d,  for  gages  7  and  8  to  avoid  a  negative  water  depth  connotation  in  the 
swash  zone,  where  d  <  0.  Rubber  lids  were  placed  on  the  14  gage  wells  that  were 
not  in  use  for  the  free  surface  measurements  in  individual  runs.  Rubber  lids  were 
also  attached  to  wave  gages  7  and  8  to  seal  the  gage  wells  used  for  the  free  surface 
measurements.  The  rubber  lids  were  used  to  provide  a  continuous,  smooth  slope 
over  the  wells  to  eliminate  disturbances  in  the  sensitive  shallow  water  and  swash 
regions.  Wave  gages  7  and  8  with  the  attached  rubber  lids,  which  were  calibrated 
before  and  after  each  test,  were  deployed  in  the  region  of  x  >  10  m  and  d  <  12.5  cm 
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Table  3.3:  Depth  of  Water  at  each  Gage  Position  for  Free  Surface  Test 


Run 

Water  Depth  d  (cm) 

Bottom  Elevation  Zh(cm) 

Gage  ^ 

Gage  5 

Gage  6 

Gage  7 

Gage  8 

3-1 

66 

65 

64 

-12.5 

-10.0 

3-2 

63 

62 

61 

-10.0 

-7.5 

msm 

59 

58 

-7.5 

57 

56 

55 

-5.0 

-2.5 

54 

53 

52 

-2.5 

0.0 

CO 

51 

50 

49 

0.0 

2.5 

48 

47 

46 

2.5 

5.0 

45 

44 

43 

5.0 

7.5 

3-9 

42 

41 

40 

7.5 

10.0 

3-10 

39 

38 

37 

10.0 

12.5 

3-11 

36 

35 

34 

12.5 

15.0 

3-12 

33 

32 

31 

15.0 

17.5 

3-13 

30 

29 

28 

17.5 

20.0 

3-14 

27 

26 

25 

20.0 

22.5 

3-15 

24 

23 

22 

20.0 

22.5 

3-16 

21 

20 

19 

22.5 

25.0 

3-17 

18 

17 

16 

22.5 

25.0 

up  to  the  upper  limit  of  wave  runup  for  each  test.  The  free  surface  measurements  at 
the  gage  wells  were  made  more  than  once  because  the  17  runs  for  each  test  allowed 
34  free  surface  measurements  using  wave  gages  7  and  8. 

3.3  Velocity  Tests 

The  velocity  measurements  for  each  of  the  three  tests  were  made  at  mid¬ 
depth,  z  =  (—d/2),  at  the  horizontal  locations  of  d  =  10,  12.5  and  (15  +  3j)  cm 
with  j  =  1,  2,  . . . ,  17  where  2  =  vertical  coordinate  taken  to  be  positive  upward 
with  2  =  0  at  SWL.  The  velocity  test  runs  are  numbered  18  through  34  for  each  test 
set.  The  fixed  wave  gages  1-3  were  also  deployed  to  ensure  the  repeatability  of  the 
wave  conditions  during  the  velocity  measurements  for  each  test.  The  free  surface 
elevations  at  the  locations  of  the  velocity  measurements  were  measured  separately 


in  the  17  runs  described  above.  A  3D  side-looking  probe  attached  to  an  acoustic 
doppler  velocimeter  was  used  to  measure  the  three  velocity  components  in  the  region 
d  =  18-66  cm,  whereas  a  2D  side-looking  probe  was  employed  to  measure  the  two 
horizontal  velocity  components  in  the  region  d  =  10-12.5  cm.  The  velocity  mea¬ 
surements  for  test  1  were  performed  twice  for  the  range  of  d  =  18-66  cm,  with  the 
2D  probe  used  from  d  =  18-27  cm  in  the  second  test  set,  to  check  the  repeatability 
of  the  velocity  measurements  and  the  sensitivity  of  the  3D  probe  in  shallow  water. 
The  cross-shore  horizontal  and  vertical  velocity  components  were  analyzed  in  the 
same  way  as  the  free  surface  elevations.  The  measured  cross-tank  horizontal  veloc¬ 
ity  components  were  negligible.  The  measured  vertical  velocity  components  were 
small  relative  to  the  corresponding  horizontal  velocity  components  and  affected  more 
by  the  turbulence  generated  by  wave  breaking.  Consequently,  only  the  horizontal 
velocity  components  are  presented  herein. 

3.4  Repeatability  of  Tests 

After  completing  the  17  runs  for  the  free  surface  and  velocity  measurements 
for  the  three  test  sets,  the  repeatability  of  the  incident  and  reflected  waves  was 
checked  for  each  set.  The  separated  incident  and  reflected  wave  spectra  at  the 
location  of  wave  gage  1  were  plotted  together  for  the  17  runs  of  each  free  surface 
and  velocity  test.  The  plotted  17  spectra  for  each  individual  test  set  were  almost 
identical,  as  shown  in  Figure  3.2  -  3.7. 

Comparisons  were  also  made  of  the  spectral  parameters  such  as  the  incident 
wave  spectral  peak  period,  Tp,  the  spectral  estimate  of  the  incident  significant  wave 
height,  Hmo  =  4^/mo,  with  mo  =  zero-moment  of  the  incident  wave  spectrum,  and 
the  average  reflection  coefficient,  R  =  y(jn^r I with  (mo)r-  =  zero-moment  of 
the  reflected  wave  spectrum.  In  addition,  the  separated  incident  wave  time  series 
were  analyzed  using  a  zero-upcrossing  method  to  obtain  the  time  series  parameters 
such  as  the  significant  wave  period  and  the  significant  wave  height  Hs-  Tables  3.4 
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Figure  3.2:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Free 
Surface  Test  1 


Figure  3.3:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Ve¬ 
locity  Test  1 


and  3.5  list  the  average  values  of  Tp,  Hmo  and  R  for  the  17  runs  for  each  of 

the  three  free  surface  tests.  The  deviations  from  the  average  values  of  T*,  Tp 
and  Hmo  among  the  17  runs  were  generally  less  than  1%  and  did  not  exceed  4%. 
The  reflection  coefficient  R  varied  more  and,  for  example,  R  was  in  the  range  R  = 
0.14-0.15  for  test  1,R  =  0.13-0.19  for  test  2,  and  R  =  0.16-0.17  for  test  3  as  shown 
in  Table  3.4.  Finally,  the  surf  similarity  parameter  ^  (Battjes  1974)  based  on  Tp 
and  Hmo  for  the  1:16  slope  is  listed  in  Tables  3.4  and  3.5  to  indicate  the  transition 
from  mostly  spilling  breakers  for  test  1  to  predominantly  plunging  breakers  for  test 
3.  Tables  3.6  -  3.11  show  the  values  of  T„  H^,  Hmo,  Hrms  and  R  for  each  run  in 
test  set  1  through  test  set  3  for  the  free  surface  and  velocity  tests. 

Additionally,  the  values  for  wave  setup  rj,  root-mean-square  wave  height  Hrms, 
skewness,  and  kurtosis  were  compared  for  gages  1-3  for  each  run  in  each  test  set. 
The  values,  shown  in  Tables  3.12  -  3.17,  are  consistent  through  each  test  set  and 
were  generally  within  small  percentage  points  of  error. 


Figure  3.4:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Free 
Surface  Test  2 


Figure  3.6:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Free 
Surface  Test  3 


Figure  3.5:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Ve¬ 
locity  Test  2 


Figure  3.7:  Smoothed  Power  Density 
Spectrum;  Incident  and 
Reflected  Waves  for  Ve¬ 
locity  Test  3 
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Ikewness,  and  Kurtosis  for  Velocity  Test  1 
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Chapter  4 


FREE  SURFACE  ELEVATIONS 


The  wave  setup  the  standard  derivation  cr,  the  skewness  s,  and  the  kurtosis 
K  for  each  of  the  measured  time  series  of  the  free  surface  elevation,  t],  of  the  325  s 
duration  sampled  at  the  rate  of  20  Hz  are  computed  under  the  assumption  of  equiv¬ 
alency  of  probabilistic  and  time  averaging.  These  time  series  include  the  effects  of 
reflected  waves.  The  root-mean-square  wave  height  H^ms  is  obtained  using  equation 
(2.8).  The  normalized  free  surface  elevation,  =  {tj  —  tJ)  j a,  is  then  calculated  us¬ 
ing  the  measured  values  of  rj  and  a.  For  each  of  the  time  series  of  77,,  the  probability 
density  function,  /(?/.),  is  computed  using  20  bins  in  the  range  determined  by  the 
measured  maximum  and  minimum  values  of  7;*. 

4.1  Probability  Density  Functions  for  Gages  1-3  for  Tests  1-3 

Figure  4.1  shows  the  measured  probability  density  functions  at  wave  gages 
1-3  for  the  17  runs  in  free  surface  tests  1-3.  The  exponential  gamma  distribution, 
fiv*)^  given  by  equation  (2.1)  is  also  plotted  in  Figure  4.1  for  the  measured  minimum 
and  maximum  values  of  s  for  each  test.  The  ordinate  for  the  figures  plotted  for  /(t/*) 
is  the  normalized  free  surface  elevation  77*  varying  vertically.  Figure  4.1  is  another 
indicator  of  the  degree  of  repeatability  for  the  17  runs  in  each  test,  as  the  different 
runs  for  each  test  set  axe  virtually  indistinguishable.  It  is  noted  that  separate  figures 
for  each  wave  gage  indicate  slightly  better  repeatability.  The  exponential  gamma 
distribution  changes  little  for  the  measured  range  of  s  for  each  test  and  is  in  good 


Wave  Gages  1,  2,  and  3 


Figure  4.1:  Measured  and  Computed  Probability  Distributions  for  17  Runs  at 
Wave  Gages  1-3  for  Test  Sets  1-3 

agreement  with  the  measured  distributions  except  for  the  scatter  of  data  points  at 
the  peak  of  /(t?*). 

4.2  Spacial  Variations  of  Statistical  Values  in  Free  Surface  Elevations 

The  measured  cross-shore  variations  of  rj,  and  s  for  tests  1-3  are 

shown  in  Figures  4.2  -  4.4,  where  the  horizontal  coordinate  a:  =  0  at  the  still  water 
depth  d  =  75  cm.  The  mean  water  depth  h  is  given  hy  h  =  {rj  —  Zb)  where 
the  bottom  elevation  Zb  corresponding  to  the  1:16  slope  is  indicated  by  the  solid 
line  in  these  figures  to  show  that  the  wave  setup  rf  becomes  tangential  to  the  1:16 
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Test  1 


x(m) 

Figure  4.2:  Measured  Cross-Shore  Variations  of  H^jas/h,  and  s  for  Test  1 

slope.  The  present  data  analysis  is  limited  to  the  region  of  h  >  0.4  cm  because 
the  statistical  analysis  becomes  less  reliable  in  the  region  of  extremely  small  depth 
affected  by  infrequent  wave  uprush  and  downrush.  The  values  of  rj,  H^mslh  and 
s  at  the  locations  of  wave  gages  1-3  are  plotted  for  all  17  runs  in  each  test.  The 
data  from  the  repeated  measurements  at  the  locations  of  the  gage  wells  are  plotted 
individually  without  averaging  to  indicate  that  the  visible  scatter  of  data  points  is 
limited  to  the  swash  zone. 

The  measured  values  of  i/rms/h  in  Figures  4.2  -  4.4  show  a  gradual  increase 
and  then  rise  rapidly  near  the  still  water  shoreline.  The  gradual  landward  increase 
of  Hrmsih  is  mostly  caused  by  the  landward  decrease  of  h  where  the  measured  values 
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Test  2 


x(m) 

Figure  4.3:  Measured  Cross-Shore  Variations  of  Hrmsih,  and  s  for  Test  2 

of  Hrms  are  approximately  constant  in  the  shoaling  zone  and  decrease  gradually  in 
the  surf  and  swash  zones.  The  measured  values  of  H^/h  on  the  1:16  slope  do 
not  approach  constant  in  the  inner  surf  zone,  unlike  the  field  data  on  gently  sloping 
concave  beaches  obtained  by  Thornton  and  Guza  (1982,  1983).  The  measured  values 
of  Hrmsih  approach  the  intuitive  upper  limit  of  =  \/8  in  the  swash  zone 

predicted  by  equation  (2.7)  with  (2.8),  although,  the  exact  upper  limit  is  difficult  to 
pinpoint  because  of  the  rapid  increase  and  scatter  of  Hrmsih  as  h  approaches  zero. 
This  implies  that  the  landward  decrease  of  ifrms  is  more  gradual  than  that  of  h  in 
the  swash  zone. 
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Test  3 


x(m) 

Figure  4,4:  Measured  Cross-Shore  Variations  of  s  for  Test  3 

Figures  4.2  -  4.4  show  that  the  skewness  6  increases  landward  and  then  de- 
creases  toward  the  still  water  shoreline  as  expected  from  available  data  [e.g.,  Mase 
and  Kobayashi  (1991);  Raubenheimer  et  al.  (1995)].  However,  the  skewness  in¬ 
creases  again  in  the  swash  zone  and  approaches  the  upper  limit  of  s  =  2  correspond¬ 
ing  to  the  exponential  distribution  given  by  equation  (2.6).  This  upper  limit  of  s  =  2 
may  not  be  very  accurate  but  implies  that  the  probability  distribution  of  the  free 
surface  elevation  is  strongly  affected  by  the  bottom  elevation.  It  is  also  noted  that 
the  detailed  cross-shore  variations  of  s  are  noticeably  different  for  the  three  tests 
whose  incident  wave  characteristics  are  listed  in  Table  3.4.  Complete  listings  of  the 
data  points  plotted  in  these  figures  can  be  found  in  Tables  A.l  -  A.3  in  Appendix 
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A. 


4.3  Measured  and  Theoretical  Probability  Density  Functions 

The  measured  probability  density  function  of  ?/*  in  the  region  of  d  <  66  cm 
was  compared  individually  with  the  exponential  gamma  distribution  given  by  equa¬ 
tion  (2.1)  together  with  the  measured  value  of  s.  The  agreement  was  good  in  the 
shoaling  and  surf  zones,  and  fair  in  the  swash  zone.  Figures  4.5  ~  4.7  show  the 
degree  of  the  overall  agreement  for  tests  1—3,  respectively.  The  region  of  d  <  66  cm 


Test  1 


Figure  4.5:  Measured  and  Computed  Probability  Distributions  in  Four  Zones  for 
Test  1 


is  separated  into  four  zones  for  each  test  on  the  basis  of  the  cross-shore  variations 
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Figure  4.6:  Measured  and  Computed  Probability  Distributions  in  Four  Zones  for 
Test  2 


of  s  shown  in  Figures  4.2  -  4.4.  The  four  zones  roughly  correspond  to  the  landward 
increase  zone  of  little  wave  breaking,  the  peak  zone  of  intense  wave  breaking,  the 
trough  zone  seaward  of  the  still  water  shoreline,  and  the  rapid  increase  zone  above 
the  still  water  shoreline.  Figures  4.5  -  4.7  indicate  the  range  of  the  still  water  depth 
d  or  the  bottom  elevation  Zb  above  SWL  for  each  zone  where  d  =  —Zb-  The  measured 
minimum  and  maximum  values  of  s  in  each  zone  are  indicated  in  these  figures.  The 
exponential  gamma  distributions  based  on  these  minimum  and  maximum  values  of  s 
are  plotted  for  each  zone  so  that  the  data  points  would  fall  between  the  two  theoret¬ 
ical  curves  if  the  agreement  were  perfect.  In  the  swash  zone,  the  exponential  gamma 
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Test  3 


Figure  4.7:  Measured  and  Computed  Probability  Distributions  in  Four  Zones  for 
Test  3 


distribution  is  limited  by  the  exponential  distribution  with  s  =  2  given  by  equation 
(2.6).  In  the  rapid  increase  zone  above  the  still  water  shoreline,  the  agreement  is 
better  for  the  large  values  of  t/,  corresponding  to  the  large  free  surface  elevations. 
This  is  partialy  because  the  theoretical  distributions  for  the  large  values  of  are 
not  very  sensitive  to  the  skewness,  s.  In  addition,  the  free  surface  measurements 
landward  of  the  still  water  shoreline  are  less  accurate  when  a  thin  layer  of  water 
remains  on  the  1:16  slope.  Separate  figures  for  each  wave  gage  are  plotted  in  Figures 
4.9  -  4.55,  and  indicate  even  better  agreement.  In  these  figures,  each  gage  position 
is  plotted  separately.  In  the  area  where  d  <  12.5  cm,  repeated  tests  are  plotted 
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together. 

Finally,  the  measured  values  of  the  kurtosis  K  are  compared  with  the  theoret¬ 
ical  values  of  K  based  on  equation  (2.4)  with  (2.3)  using  the  corresponding  measured 
values  of  s.  Figure  4.8  shows  all  the  measured  values  of  {K  —  3)  as  a  function  of  s 
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Figure  4.8;  Measured  and  Computed  Kurtosis,  Ji ,  as  a  Function  of  Measured 
Skewness,  s 


Test  1 , 2,  and  3 


in  comparison  to  the  theoretical  curve  shown  in  Figure  2.1  where  (K  —  3)  =  6  for 
5  >  2  is  assumed.  The  relation  between  K  and  s  based  on  the  exponential  gamma 
distribution  tends  to  overpredict  K  somewhat  for  s  <  2  and  can  not  explain  the 
increase  of  K  for  s  >  2. 


Figure  4.9:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  66,  65,  64,  and  63cm  for  Test  1 
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Figure  4.10:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  62,  61,  60,  and  59cm  for  Test  1 
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Figure  4.12:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  54,  53,  52,  and  51cm  for  Test  1 
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Figure  4.13:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  50,  49,  48,  and  47cm  for  Test  1 


Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  —  42,  41,  40,  and  39cm  for  Test  1 
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Figure  4.18:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  30,  29,  28,  and  27cm  for  Test  1 
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Figure  4.24:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  66,  65,  64,  and  63cm  for  Test  2 
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Figure  4.27:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  54,  53,  52,  and  51cm  for  Test  2 
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Figure  4.28:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  —  50,  49,  48,  and  47cm  for  Test  2 
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Figure  4.30:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  42,  41,  40,  and  39cm  for  Test  2 
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Figure  4.31:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  38,  37,  36,  and  35cm  for  Test  2 
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Figure  4.34:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  26,  25,  24,  and  23cm  for  Test  2 
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Figure  4.40:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  66,  65,  64,  and  63cm  for  Test  3 


Figure  4.43: 


Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  54,  53,  52,  and  51cm  for  Test  3 
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Figure  4.45:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  46,  45,  44,  and  43cm  for  Test  3 
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Figure  4.52:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  18,  17,  and  16cm  for  Test  3 


Figure  4.53:  Measured  and  Computed  Probability  Distributions  at  Water  Depth 
d  =  12.5,  10,  7.5,  and  5cm  for  Test  3 
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Figure  4.55:  Measured  and  Computed  Probability  Distributions  at  Bottom  Ele 
vation  =  7.5,  10,  12.5,  and  15cm  for  Test  3 


Chapter  5 


CROSS  SHORE  HORIZONTAL  VELOCITIES 

The  mean  velocity  ¥,  the  standard  deviation  the  skewness  s,,,  and  the  kur- 
tosis  Ku  for  each  of  the  time  series  of  the  cross-shore  horizontal  velocity  u  measured 
at  the  mid-depth  below  SWL  are  computed  under  the  assumption  of  equivalency 
of  the  probabilistic  and  time  averaging.  The  horizontal  velocity  u  is  taken  to  be 
positive  landward  and  the  undertow  velocity  u  is  negative  herein.  The  normalized 
horizontal  velocity,  u.  =  (u  -  u)  fau,  is  then  calculated  using  the  measured  values  of 
u  and  cr„.  The  probability  density  function,  /(w*),  for  each  of  the  time  series  of  u,  is 
obtained  and  compared  with  the  exponential  gamma  distribution  given  by  equation 
(2.1)  with  (2.2)  and  (2.3)  where  r/,  and  s  are  replaced  by  u,  and  respectively. 
The  skewness  and  kurtosis  are  the  same  for  u  and  u,. 

5.1  Horizontal  Velocity  and  Linear  Theory 

To  obtain  approximate  relationships  between  the  free  surface  and  horizontal 
velocity  statistics,  local  nonlinearity  may  be  neglected  and  linear  progressive  long¬ 
wave  theory  may  be  assumed  to  be  approximately  valid  locally  even  inside  surf  zones 
(Guza  and  Thornton  1980).  The  oscillatory  components  of  the  horizontal  velocity 
and  free  surface  elevation  may  then  be  related  in  the  following  form: 


in  which  g  =  gravitational  acceleration  and  h  =  mean  water  depth  given  hy  h  = 
{d  +  If)  with  d  =  still  water  depth.  If  equation  (5.1)  holds,  the  following  relationships 
can  be  derived 


in  which  cr,  s  and  K  are  the  standard  deviation,  skewness,  and  kurtosis  of  the 
free  surface  elevation  rj.  Furthermore,  the  probability  density  functions  of  m*  = 
(u  —  u)  jay,  and  =  {v  ~v)  I ^  become  approximately  the  same  because  u*  ~ 
under  the  assumption  of  (5.1)  with  cr„  ~  yjgjha. 

In  order  to  estimate  the  undertow  velocity  n,  the  horizontal  velocity  u  at  the 
mid-depth  is  assumed  to  be  approximately  equal  to  the  depth- averaged  velocity  U. 
The  time-averaged  continuity  equation,  hU  =  0,  for  the  impermeable  beach  yields 
(Kobayashi  et  al.  1989) 

tr  =-(!)■' (,-?)(£/- £7)  (5.3) 

in  which  use  is  made  of  =  (d  -f  g).  Substitution  of  u  ~  1/  and  equation  (5.1)  into 
(5.3)  yields 


5.2  Spacial  Variations  of  Statistical  Values  in  Horizontal  Velocity 

Figures  5.1  -  5.3  show  the  cross-shore  variations  of  the  measured  values  of  u, 
(T„,  Sy  and  Ky  in  the  still  water  depth  d  =  10-66  cm  for  tests  1-3,  with  the  values 
also  listed  in  Tables  5.1  -  5.3.  The  two  data  points  at  each  cross-shore  location, 
except  for  the  landward  two  locations  in  Figure  5.1,  correspond  to  the  repeated 
velocity  measurements  for  test  1.  The  repeatability  of  the  velocity  measurements  is 
good  except  for  the  undertow  velocity  u  in  the  shallow  water.  The  computed  values 
of  u,  (Xy,  Sy,  and  Ky  at  each  cross-shore  location  shown  in  Figures  5.1  -  5.3  are  based 
on  equations  (5.2)  and  (5.4)  together  with  the  measured  values  oi  h  =  [d  +  rj),  a,  s, 


and  K  at  the  same  location. 


Figure  5.1:  Measured  and  Computed  Cross-Shore  Variations  of  u,  <Tu,  and  A„ 
for  Test  1 

The  fair  (within  about  20%  errors)  agreement  between  the  measured  and 
computed  standard  deviation  <Tu  in  Figures  5.1  -  5.3  is  expected  from  the  field  data 
intercomparisons  made  by  Guza  and  Thornton  (1980).  For  d  =  10-66  cm  with  Lp 
=  linear  wavelength  based  on  Tp,  Lpjd  =  5-14,  10-27,  and  18—46  for  tests  1,  2, 
and  3,  respectively.  The  long-wave  assumption  may  not  be  appropriate  for  test 
1.  It  is  noted  that  Guza  and  Thornton  (1980)  limited  their  analysis  to  the  wind 
wave  frequencies,  whereas  the  measured  time  series  without  any  filtering  are  used 
herein  because  it  is  difficult  to  separate  wind  wave  and  low-frequency  components 
for  practical  applications.  The  fair  agreement  between  the  measured  and  computed 
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Figure  5.2:  Measured  and  Computed  Cross- Shore  Variations  of  u,  (7„,  and  AV 
for  Test  2 


undertow  u  is  unexpected  because  equation  (5.4)  based  on  (5.1)  does  not  account  for 
the  additional  water  volume  flux  due  to  a  roller  which  increases  the  magnitude  of  u  as 
shown  by  Svendsen  (1984)  for  regular  waves.  The  underprediction  of  the  magnitude 
of  u  in  shallow  water  for  tests  2  and  3  may  be  explained  by  this  roller  effect.  The 
landward  increases  of  \Ti\  and  in  Figures  5.1  -  5.3  are  qualitatively  consistent  with 
those  based  on  the  depth- averaged  velocity  U  computed  by  Kobayashi  and  Karjadi 
(1996)  using  the  finite-amplitude  shallow- water  equations.  In  view  of  equations 
(5.2)  and  (5.4),  the  landward  increases  of  |u|  and  cr„  may  simply  be  explained  by 
the  landward  increase  of  H„ns/h  shown  in  Figures  4.2  -  4.4.  In  conclusion,  the 
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X  (m) 

Figure  5.3:  Measured  and  Computed  Cross-Shore  Variations  of  u,  <Ju,  Su,  and  Ku 
for  Test  3 


intercomparisons  by  Guza  and  Thornton  (1980)  and  Figures  5.1  -  5.3  indicate  that 
the  local  application  of  linear  long-wave  theory  is  approximately  valid  for  the  wave 
quantities  up  to  the  second  order  or  moment  in  the  shoaling  and  surf  zones. 

On  the  other  hand,  Figures  5.1  -  5.3  show  that  the  skewness  and  kurtosis 
representing  the  third  and  fourth  moments  are  not  predicted  well  by  linear  long- wave 
theory.  The  skewness  Su  of  the  horizontal  velocity  u  at  the  mid-depth  is  always  less 
than  the  skewness  s  of  the  free  surface  elevation  ?/.  This  implies  that  nonlinearity 
decreases  downward  from  the  free  surface.  The  difference  (5  —  s„)  tends  to  increase 
landward  but  fluctuates,  partly  because  of  the  difficulty  in  measuring  velocities  in 
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very  shallow  water.  On  the  average,  the  values  of  (s  —  5u)  are  on  the  order  of  0.5. 
The  measured  values  of  Su  for  the  three  tests  exhibit  noticeable  differences.  The 
skewness  in  the  shoaling  regions  in  Figures  5.1  -  5.3  increases  with  the  increase 
of  the  surf  similarity  parameter  (  where  ^  =  0.29,  0.47  and  0.72  for  tests  1,  2  and 
3,  respectively,  as  listed  in  Table  3.5.  The  skewness  Su  in  this  shoaling  region  for 
test  1  is  small  and  negative,  which  implies  that  the  cross-shore  velocities  are  skewed 
slightly  seaward.  The  kurtosis,  Ku,  of  the  mid-depth  horizontal  velocity,  u,  and  the 
kurtosis,  IF,  of  the  free  surface  elevation,  t;,  exhibit  trends  similar  to  and  s  except 
that  the  measured  kurtosis,  Ku,  in  shallow  water  fluctuates  widely  and  may  not  be 
very  reliable.  It  is  noted  that  the  measured  values  of  Su  and  Ku  for  these  three  tests 
are  qualitatively  similar  to  the  field  data  of  Guza  and  Thornton  (1985). 

5.3  Measured  and  Theoretical  Probability  Density  Functions 

The  measured  probability  density  function,  /(u*),  is  compared  individually 
with  the  exponential  gamma  distribution  based  on  the  measured  skewness,  The 
agreement  is  good  except  for  slight  deviations  at  the  peaks  of  /(u*).  To  show 
the  degree  of  the  overall  agreement  for  tests  1-3,  the  measured  distributions  are 
separated  into  two  groups  on  the  basis  of  the  measured  skewness  Figure  5.4 
shows  the  measured  distributions  with  Su  =  (-0.25)  -  0.36  in  comparison  with  the 
exponential  gamma  distributions  with  =  0  and  0.36  where  >  0  is  assumed  to  be 
consistent  with  the  assumption  of  s  >  0  made  in  (2.1).  The  measured  distributions 
with  Su  —  0.39  -  0.67  are  compared  with  the  exponential  gamma  distributions  with 
Su  =  0.39  and  0.67.  The  exponential  gamma  distributions  do  not  deviate  much 
within  the  specified  range  of  Su-  Figure  5.4  indicates  that  most  of  the  data  points 
fall  between  the  two  theoretical  curves,  apart  from  the  scatter  near  the  peak  of 
/(u,).  The  agreement  is  similar  to  those  shown  in  Figures  4.5  -  4.7  for  the  free 
surface  elevations.  It  is  not  certain  whether  the  exponential  gamma  distribution  is 


applicable  to  horizoBtal  velocities  near  the  still  water  shoreline  for  lack  of  data.  The 
individual  comparisons  are  shown  in  Figures  5.5  -  5.19  for  completeness. 
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■0.125  13.31  3.89  -12.23  -9.76  29.95  33.40  0.35  1.02  3.85  4.35 

-0.10  11.31  3.47  -12.00  -9.92  29.37  32.32  0.26  0.80  3.86  3.61 


Figure  5.5:  Measured  and  Computed  Probability  Distributions  for  Horizontal  Ve 
locity  at  Water  Depth  d  =  66,  63,  60,  and  57cm  for  Test  1 


Test  1 


Figure  5.9:  Measured  and  Computed  Probability  Distributions  for  Horizontal  Ve 
locity  at  Water  Depth  d  =  18,  12.5,  and  10cm  for  Test  1 


Figure  5.11:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  54,  51,  48,  and  45cm  for  Test  2 
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Test  2 


Figure  5.13:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  —  30,  27,  24,  and  21cm  for  Test  2 


104 


Test  2 


Figure  5.14:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  18,  12.5,  and  10cm  for  Test  2 


Tests 


Figure  5.15:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  66,  63,  60,  and  57cm  for  Test  3 
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3 


Figure  5.16:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  54,  51,  48,  and  45cm  for  Test  3 
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Tests 


Figure  5.17:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  42,  39,  36,  and  33cm  for  Test  3 


Tests 


Figure  5.18:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  30,  27,  24,  and  21cm  for  Test  3 


Tests 


Figure  5.19:  Measured  and  Computed  Probability  Distributions  for  Horizontal 
Velocity  at  Water  Depth  d  =  18,  12.5,  and  10cm  for  Test  3 


Chapter  6 


CONCLUSION 

Three  irregular  wave  tests  were  conducted  on  a  1:16  smooth  slope  to  investi¬ 
gate  the  detailed  cross-shore  variations  of  the  probability  distributions  and  statistics 
of  the  free  surface  elevations  and  mid-depth  horizontal  velocities  in  the  shoaling,  surf 
and  swash  zones.  The  exponential  gamma  distribution  is  adopted  in  this  study  be¬ 
cause  it  becomes  Gaussian  for  zero  skewness  and  exponential  for  the  skewness  of 
two.  The  exponential  gamma  distribution  with  the  measured  mean,  standard  devia¬ 
tion,  and  skewness  is  shown  to  be  capable  of  describing  all  the  measured  probability 
distributions  in  a  unified  manner. 

The  lower  limit  of  the  exponential  distribution  of  the  free  surface  elevation, 
T],  is  assumed  to  be  imposed  by  the  bottom  elevation  in  the  swash  zone.  This 
assumption  yields  =  t/S,  where  h  =  mean  water  depth,  and  H^ms  =  root- 

mean-square  wave  height  defined  as  i/rms  —  with  a  =  standard  deviation  of 
1].  The  measured  values  of  gradually  increase  in  the  surf  zone  and  rapidly 

approach  this  upper  limit  of  approximately  three  in  the  swash  zone.  The  measured 
skewness,  s,  of  77  increases  landward  and  then  decreases  toward  the  still  water  shore¬ 
line  but  increases  again  in  the  swash  zone  to  rapidly  approach  the  upper  limit  of 
s  =  2.  The  measured  kurtosis,  K,  of  77  can  also  be  explained  by  the  relation  between 
K  and  s  derived  from  the  exponential  gamma  distribution.  The  major  conclusion 
is  that  the  probability  distribution  and  statistics  of  77  in  the  swash  zone  are  very 
different  from  those  in  the  surf  zone.  The  simple  relationships  of  Hims  —  VS  h  and 
5  =  2  are  useful  in  predicting  the  wave  motion  in  the  swash  zone. 


Ill 


The  mean  (undertow)  u,  and  the  standard  deviation,  (7„,  of  the  horizontal 
velocity,  u,  measured  at  the  mid-depth  below  SWL  in  the  shoaling  and  surf  zones  are 
expressed  in  terms  of  h  and  Hrms  using  linear  progressive  long-wave  theory.  These 
simple  relationships  together  with  the  measured  values  of  h  and  Hms  can  explain 
the  measured  cross-shore  variations  of  u  and  ciu  fairly  accurately.  The  skewness, 
Su,  of  u  is  found  to  be  always  smaller  than  the  corresponding  skewness,  s,  of  i], 
indicating  the  vertical  decrease  of  nonlinearity.  This  can  not  be  explained  by  linear 
long- wave  theory.  The  kurtosis,  Ku,  oi  u  exhibits  a  similar  trend  but  the  measured 
values  are  less  reliable  in  shallow  water. 

To  utilize  the  exponential  gamma  distribution  and  the  simple  relationships 
presented  herein  for  practical  applications,  a  time- averaged  model  will  need  to  be 
developed  to  predict  the  cross-shore  variations  of  h,  a,  s,  and  s„.  The  equivalency 
of  probabilistic  and  time  averaging  will  allow  one  to  estimate  various  even  and  odd 
moments  associated  with  time-averaged  cross-shore  sediment  transport  rates  as  was 
done  by  Guza  and  Thornton  (1985)  using  the  Gaussian  distribution  which  predicts 
zero  odd  moments. 
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Appendix  A 


COMPLETE  FREE  SURFACE  STATISTICAL  DATA  FOR 

TEST  SETS  1-3 


Complete  free  surface  statistical  data  is  given  below  for  each  gage  location 
for  test  1.  Each  gage  is  indexed  by  its  x  position  and  the  corresponding  water  depth 
where  =  —d. 

Table  A.l;  Test  1  Gage  Data 


X 

n 

h 

Hrms 

s 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

datcL 

theory 

12.44 

0.36 

3.23 

3.26 

0.38 

-72.6 

12.49 

0.38 

3.20 

3.29 

0.50 

0.36 

3.17 

3.25 

1.44 

warn 

0.24 

3.02 

3.11 

1.60 

-65.0 

0.24 

2.99 

3.11 

1.76 

-64.0 

0.22 

2.99 

3.10 

1.92 

-63.0 

-0.07 

62.93 

0.28 

3.01 

3.16 

2.08 

-62.0 

-0.08 

61.92 

0.31 

3.07 

3.19 

2.24 

-61.0 

— 

0.31 

3.01 

3.19 

2.40 

-60.0 

0.33 

3.02 

3.22 

-59.0 

0.35 

3.04 

mmm 

0.36 

3.05 

Table  A.l:  (Test  1  Gage  Data,  continued) 


X 

Zb 

V 

h 

Hfins 

s 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

-57.0 

-0.28 

56.72 

11.61 

0.37 

3.07 

3.27 

-56.0 

-0.07 

55.93 

11.62 

0.37 

3.08 

3.27 

-55.0 

-0.09 

54.91 

11.63 

0.37 

3.10 

3.27 

-54.0 

-0.05 

53.95 

11.74 

0.42 

3.19 

3.35 

3.52 

-53.0 

-0.08 

52.92 

11.67 

0.40 

3.15 

3.33 

3.68 

-52.0 

-0.10 

51.90 

11.56 

0.39 

3.14 

3.30 

3.84 

-51.0 

-0.13 

50.87 

11.55 

0.44 

3.18 

3.39 

4.00 

-50.0 

-0.06 

49.94 

11.43 

0.43 

3.15 

3.36 

4.16 

-49.0 

-0.15 

48.85 

11.32 

0.43 

3.16 

3.37 

4.32 

-48.0 

47.99 

11.43 

0.46 

3.15 

3.42 

4.48 

-47.0 

46.89 

11.49 

0.46 

3.14 

3.42 

4.64 

-46.0 

-0.09 

45.91 

11.55 

0.44 

3.13 

3.38 

4.80 

-45.0 

-0.00 

45.00 

11.75 

0.49 

3.23 

3.48 

4.96 

-44.0 

-0.10 

43.90 

11.68 

0.49 

3.21 

I^QII 

5.12 

-43.0 

-0.10 

42.90 

11.53 

0.49 

3.23 

3.48 

5.28 

-42.0 

-0.07 

41.93 

11.42 

0.53 

3.27 

3.55 

5.44 

-41.0 

-0.14 

40.86 

11.23 

0.54 

3.26 

3.58 

5.60 

-40.0 

39.85 

11.02 

0.52 

3.18 

3.53 

5.76 

-39.0 

39.03 

11.07 

0.56 

3.24 

3.61 

5.92 

-38.0 

37.88 

11.02 

0.57 

3.22 

3.64 

6.08 

-37.0 

-0.18 

36.82 

10.98 

0.56 

3.20 

3.61 

6.24 

-36.0 

-0.20 

35.80 

10.98 

0.60 

3.27 

3.72 

6.40 

-35.0 

-0.24 

34.76 

10.79 

0.60 

3.23 

3.72 
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Table  A.l;  (Test  1  Gage  Data,  continued) 


X 

V 

h 

Hrms 

5 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

6.56 

-34.0 

-0.26 

33.74 

10.86 

0.62 

3.22 

3.77 

6.72 

-33.0 

-0.10 

32.90 

10.94 

0.69 

3.33 

3.95 

6.88 

-32.0 

-0.15 

31.85 

10.77 

0.72 

3.34 

4.01 

7.04 

-31.0 

-0.15 

30.85 

10.79 

0.74 

3.42 

4.07 

7.20 

-30.0 

-0.09 

29.91 

10.85 

0.81 

3.54 

4.28 

7.36 

-29.0 

-0.14 

28.86 

10.61 

0.80 

3.42 

4.24 

7.52 

-28.0 

-0.13 

27.87 

10.63 

0.78 

3.31 

4.18 

-27.0 

26.95 

10.60 

0.82 

3.44 

4.32 

7.84 

-26.0 

25.84 

10.44 

0.80 

3.32 

4.25 

8.00 

-25.0 

-0.13 

24.87 

10.34 

0.77 

3.22 

4.17 

8.16 

23.93 

10.44 

0.86 

3.42 

4.44 

8.32 

22.94 

10.29 

0.87 

3.37 

4.46 

8.48 

-22.0 

-0.09 

21.91 

10.06 

0.85 

3.29 

4.40 

8.64 

-21.0 

-0.05 

20.95 

10.06 

0.93 

3.45 

4.67 

8.80 

-20.0 

0.01 

20.01 

9.86 

0.95 

3.49 

4.71 

8.96 

-19.0 

0.05 

19.05 

9.69 

0.97 

3.59 

4.80 

9.12 

-18.0 

0.06 

18.06 

9.41 

0.99 

3.63 

9.28 

^^1 

0.08 

9.44 

my 

0.12 

10.00 

0.30 

10.00 

HOI 

0.98 

3.79 

10.40 

10.57 

6.55 

0.88 

3.72 

10.40 

-10.0 

0.52 

10.52 

6.48 

0.90 

3.97 

4.55 
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Table  A.l:  (Test  1  Gage  Data,  continued) 


X 

n 

h 

Hfrns 

5 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

-10.0 

0.45 

10.45 

6.54 

0.85 

3.80 

4.38 

10.40 

-10.0 

0.55 

10.55 

6.53 

0.83 

3.71 

4.35 

10.80 

-7.5 

0.65 

8.15 

5.61 

0.67 

3.39 

3.90 

10.80 

-7.5 

0.79 

5.61 

0.68 

3.37 

3.92 

10.80 

-7.5 

0.85 

8.35 

5.76 

0.71 

3.44 

4.00 

10.80 

0.83 

8.33 

5.74 

0.69 

3.30 

3.93 

11.20 

1.14 

6.14 

4.74 

0.57 

3.09 

3.64 

11.20 

-5.0 

1.14 

6.14 

4.74 

0.54 

2.98 

3.58 

11.20 

-5.0 

1.22 

6.22 

4.87 

0.60 

3.05 

3.71 

11.20 

-5.0 

1.21 

6.21 

4.90 

0.56 

2.91 

3.62 

11.60 

my 

1.48 

3.98 

4.19 

0.36 

2.96 

3.26 

11.60 

1.48 

3.98 

4.21 

0.41 

2.96 

3.33 

11.60 

-2.5 

1.54 

4.04 

4.27 

0.41 

3.12 

3.34 

11.60 

-2.5 

1.60 

4.10 

4.31 

1 

0.42 

3.08 

3.35 

12.00 

0.0 

1.75 

1.75 

3.17 

0.64 

3.10 

3.81 

12.00 

^^1 

1.69 

1.69 

3.19 

0.67 

3.13 

3.90 

12.00 

1.83 

1.83 

3.36 

0.57 

2.90 

3.65 

12.00 

0.0 

1.83 

2.97 

3.66 

12.40 

2.5 

5.72 

7.29 

12.40 

2.5 

3.20 

0.70 

1.90 

1.59 

5.66 

7.24 

n 

10^ 

3.16 

0.66 

2.03 

1.63 

5.73 

7.41 

12.40 

HH 

3.16 

0.66 

2.03 

1.69 

6.26 

7.68 

12.80 

5.0 

5.26 

0.26 

0.98 

2.24 

9.76 

9.00 

Table  A.l;  (Test  1  Gage  Data,  continued) 


X 

(m) 

(cm) 

V 

(cm) 

h 

(cm) 

Hrms 

(cm) 

s 

Kurtosis  K 

dOftCL 

theory 

12.80 

5.0 

5.31 

0.31 

0.96 

2.24 

9.88 

9.00 

12.80 

5.0 

5.56 

0.56 

1.26 

1.97 

8.16 

8.89 

12.80 

5.0 

5.49 

0.49 

1.32 

1.98 

8.13 

8.91 

12.80 

5.0 

5.48 

0.48 

1.30 

1.95 

7.91 

8.78 

Complete  free  surface  statistical  data  is  given  below  for  each  gage  location 
for  test  2.  Each  gage  is  indexed  by  its  x  position  and  the  corresponding  water  depth 
where  Zb  =  —d. 


Table  A. 2:  Test  2  Gage  Data 


X 

^6 

1 

h 

Hrms 

s 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

0.00 

-75.0 

-0.32 

74.68 

16.85 

0.52 

3.51 

3.55 

0.38 

-72.6 

-0.24 

72.36 

16.71 

0.60 

3.58 

3.72 

0.50 

-71.9 

-0.11 

71.79 

16.53 

0.62 

3.60 

3.76 

1.44 

-66.0 

-0.20 

65.80 

16.28 

0.68 

3.76 

3.92 

1.60 

-65.0 

! 

-0.28 

64.72 

16.16 

0.67 

3.74 

3.90 

1.76 

-64.0 

-0.23 

63.77 

15.94 

0.69 

3.84 

3.95 

1.92 

-63.0 

-0.11 

62.89 

16.17 

0.74 

3.93 

4.09 

2.08 

-62.0 

-0.16 

61.84 

16.06 

0.73 

3.89 

4.04 

2.24 

-61.0 

-0.21 

60.79 

15.81 

0.74 

3.97 

4.08 

2.40 

-60.0 

-0.07 

59.93 

16.14 

0.81 

4.13 

iiom 

-59.0 

-0.12 

58.88 

16.00 

0.78 

4.06 

4.20 

Table  A. 2:  (Test  2  Gage  Data,  continued) 


X 

Zb 

h 

Hrms 

s 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

2.72 

-58.0 

-0.12 

57.88 

15.78 

0.81 

4.16 

4.28 

2.88 

-57.0 

-0.03 

56.97 

16.32 

0.88 

4.30 

4.49 

3.04 

-56.0 

-0.08 

55.92 

16.13 

0.87 

4.27 

4.47 

3.20 

^^1 

-0.08 

54.92 

15.80 

0.89 

4.30 

4.51 

3.36 

-0.04 

53.96 

16.19 

1.00 

4.58 

4.88 

3.52 

-53.0 

-0.14 

52.86 

15.85 

0.98 

4.45 

4.82 

3.68 

-0.14 

51.86 

15.41 

1.01 

4.55 

4.92 

3.84 

1^1 

-0.06 

50.94 

15.71 

1.10 

4.89 

5.24 

4.00 

-50.0 

IQQH 

49.83 

15.38 

1.10 

4.88 

5.26 

4.16 

-49.0 

48.78 

14.92 

1.14 

5.06 

5.40 

4.32 

-48.0 

-0.09 

47.91 

15.19 

1.18 

4.48 

-47.0 

-0.19 

46.81 

14.84 

1.17 

mm 

-46.0 

-0.20 

45.80 

14.42 

1.19 

mm 

4.80 

-45.0 

-0.08 

44.92 

14.52 

1.27 

5.36 

5.90 

4.96 

-44.0 

-0.13 

43.87 

14.21 

1.23 

5.22 

5.76 

5.12 

-43.0 

-0.15 

42.85 

13.84 

1.22 

5.16 

5.71 

5.28 

-42.0 

-0.09 

41.91 

14.04 

1.29 

5.34 

5.98 

5.44 

-41.0 

-0.14 

40.86 

13.87 

1.25 

5.18 

5.84 

5.60 

-40.0 

-0.12 

39.88 

13.64 

1.23 

5.12 

5.75 

5.76 

-39.0 

-0.06 

38.94 

13.99 

1.23 

5.09 

5.75 

5.92 

-38.0 

-0.11 

37.89 

13.70 

1.22 

5.05 

5.70 

6.08 

-37.0 

13.80 

1.23 

5.09 

5.74 

6.24 

-36.0 

0.08 

36.08 

14.13 

1.28 

5.19 

5.95 

Table  A. 2;  (Test  2  Gage  Data,  continued) 


Table  A. 2:  (Test  2  Gage  Data,  continued) 


X 

V 

h 

Hrms 

s 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

10.80 

-7.5 

1.82 

9.32 

9.25 

0.85 

10.80 

-7.5 

1.78 

9.28 

9.19 

0.86 

11.20 

-5.0 

2.27 

7.27 

8.34 

0.73 

11.20 

-5.0 

2.21 

7.21 

8.28 

0.87 

4.28 

4.45 

11.60 

-2.5 

2.82 

5.32 

7.64 

0.84 

3.60 

4.38 

11.60 

-2.5 

2.78 

5.28 

7.52 

0.91 

3.97 

4.59 

12.00 

0.0 

3.32 

3.32 

6.69 

1.15 

4.13 

5.43 

12.00 

0.0 

3.39 

3.39 

6.58 

1.15 

4.37 

5.42 

12.40 

2.5 

5.12 

2.62 

5.21 

1.45 

4.94 

6.64 

12.40 

2.5 

4.57 

1 

2.07 

5.25 

1.40 

4.74 

6.43 

12.80 

5.0 

6.67 

1.67 

3.66 

1.54 

5.54 

7.03 

12.80 

5.0 

6.36 

1.36 

3.56 

1.54 

5.25 

7.03 

13.20 

7.5 

8.31 

0.81 

2.33 

2.32 

9.79 

9.00 

13.20 

7.5 

8.15 

0.65 

2.23 

2.30 

9.82 

9.00 

13.60 

10.0 

11.07 

1.07 

1.78 

3.02 

14.77 

9.00 

13.60 

10.0 

10.48 

0.48 

1.96 

2.75 

12.72 

9.00 

13.60 

10.0 

10.45 

0.45 

1.95 

2.76 

13.17 

9.00 

14.00 

12.5 

13.40 

0.90 

1.44 

1.75 

8.55 

7.93 

14.00 

12.5 

13.37 

0.87 

1.43 

1.34 

7.26 

6.20 

14.00 

12.5 

13.22 

0.72 

1.23 

2.67 

13.57 

9.00 

14.00 

12.5 

12.94 

0.44 

1.06 

2.02 

9.04 

9.00 

Complete  free  surface  statistical  data  is  given  below  for  each  gage  location 


for  test  3.  Each  gage  is  indexed  by  its  x  position  and  the  corresponding  water  depth 
where  Zb  =  —d. 


Table  A. 3:  Test  3  Gage  Data 


X 

Zh 

V 

h 

Hrms 

$ 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

-1.00 

-76.2 

-0.24 

75.96 

18.40 

0.42 

2.98 

3.36 

0.38 

-72.6 

-0.30 

72.30 

18.16 

0.50 

3.22 

3.50 

0.50 

-71.9 

-0.27 

71.63 

18.30 

0.53 

3.29 

3.57 

1.44 

-66.0 

-0.21 

65.79 

18.12 

0.58 

3.46 

3.66 

1.60 

-65.0 

-0.36 

64.65 

18.18 

0.57 

3.49 

3.64 

1.76 

-64.0 

-0.39 

63.61 

18.11 

0.59 

3.56 

3.68 

1.92 

-63.0 

-0.15 

62.85 

18.13 

0.62 

3.60 

3.76 

2.08 

-62.0 

-0.26 

61.74 

18.21 

0.62 

3.64 

3.77 

2.24 

-61.0 

-0.27 

60.73 

18.10 

0.64 

3.70 

3.80 

2.40 

-60.0 

59.81 

18.09 

0.66 

3.72 

3.86 

2.56 

-59.0 

58.73 

18.14 

0.65 

3.70 

3.83 

2.72 

-58.0 

-0.27 

57.73 

18.13 

0.69 

3.86 

3.95 

2.88 

-57.0 

-0.16 

56.84 

18.10 

0.72 

3.86 

4.01 

3.04 

-56.0 

-0.31 

55.69 

18.13 

0.70 

3.84 

3.98 

3.20 

-55.0 

-0.27 

54.73 

18.00 

0.72 

3.93 

4.03 

3.36 

-54.0 

-0.18 

53.82 

18.02 

0.77 

4.07 

4.14 

3.52 

-53.0 

0.20 

53.20 

18.01 

0.75 

4.02 

4.09 

3.68 

-52.0 

-0.28 

51.72 

17.86 

0.75 

4.06 

4.10 

3.84 

-0.19 

50.81 

17.97 

0.78 

4.17 

4.17 

4.00 

-0.35 

49.65 

17.96 

0.76 

4.13 

4.13 

4.16 

-0.32 

48.68 

17.83 

0.76 

4.12 

4.13 

Table  A. 3:  (Test  3  Gage  Data,  continued) 


X 

Zh 

V 

h 

Hrms 

S 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

4.32 

-48.0 

-0.22 

47.78 

17.90 

0.81 

4.27 

4.27 

4.48 

-47.0 

-0.29 

46.71 

17.86 

0.77 

4.10 

4.16 

4.64 

-46.0 

-0.31 

45.69 

17.75 

0.76 

4.11 

4.14 

4.80 

-45.0 

-0.23 

44.77 

18.00 

0.82 

4.22 

4.29 

4.96 

-44.0 

-0.26 

43.74 

17.98 

0.78 

4.13 

4.20 

ips 

-43.0 

-0.32 

42.68 

mam 

-42.0 

-0.22 

41.78 

im^iu 

mm 

-0.33 

40.67 

— 

19 

5.60 

-40.0 

-0.33 

39.67 

■El 

lEI 

5.76 

-39.0 

-0.23 

38.77 

18.05 

0.87 

4.32 

4.45 

-38.0 

-0.28 

37.72 

18.09 

0.84 

4.21 

4.37 

II^Q 

-0.38 

36.62 

17.96 

4.19 

4.30 

-0.15 

35.85 

18.11 

4.27 

4.54 

-0.32 

34.68 

18.14 

0.87 

4.30 

4.47 

-0.38 

33.62 

17.91 

0.85 

4.38 

4.40 

6.72 

-33.0 

-0.20 

32.80 

17.73 

0.93 

4.38 

4.66 

6.88 

-32.0 

-0.29 

31.71 

17.57 

4.28 

4.56 

7.04 

-31.0 

-0.35 

30.65 

17.28 

4.41 

4.24 

7.20 

-30.0 

-0.10 

29.90 

16.95 

4.07 

4.54 

7.36 

-29.0 

-0.21 

28.79 

16.68 

0.85 

4.12 

4.40 

7.52 

-28.0 

-0.29 

27.71 

16.52 

0.68 

4.69 

3.91 

7.68 

-27.0 

-0.01 

26.99 

mam 

0.88 

4.50 

7.84 

-26.0 

-0.12 

25.88 
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Table  A. 3:  (Test  3  Gage  Data,  continued) 


X 

Zb 

n 

h 

s 

Kurtosis  K 

(m) 

(cm) 

(cm) 

(cm) 

(cm) 

data 

theory 

8.00 

-25.0 

25.18 

15.81 

0.88 

4.18 

4.50 

8.16 

-24.0 

24.23 

15.51 

0.86 

4.03 

4.42 

0.13 

23.13 

15.18 

0.74 

4.15 

4.07 

0.21 

22.21 

14.59 

0.67 

4.16 

3.88 

8.64 

-21.0 

0.45 

21.45 

14.62 

0.80 

3.79 

4.26 

8.80 

-20.0 

0.49 

20.49 

14.53 

0.70 

3.51 

3.97 

8.96 

-19.0 

0.53 

19.53 

13.79 

0.64 

3.32 

3.82 

9.12 

-18.0 

18.75 

13.48 

0.62 

3.31 

3.77 

9.28 

-17.0 

17.81 

13.31 

0.60 

3.08 

3.71 

9.44 

-16.0 

0.96 

16.96 

12.84 

0.62 

3.14 

3.75 

10.00 

-12.5 

13.93 

11.96 

0.56 

3.06 

3.61 

10.40 

-10.0 

1.85 

11.85 

11.56 

0.62 

3.20 

3.77 

-10.0 

2.01 

12.01 

11.53 

0.62 

3.20 

3.76 

-7.5 

2.25 

9.75 

10.96 

0.57 

3.10 

3.64 

10.80 

-7.5 

2.39 

9.89 

11.07 

0.57 

3.08 

3.64 

11.20 

-5.0 

2.82 

7.82 

10.52 

0.60 

3.24 

3.71 

11.20 

-5.0 

2.97 

7.97 

10.50 

0.60 

3.27 

3.71 

11.60 

-2.5 

3.58 

6.08 

9.37 

0.82 

3.37 

4.31 

11.60 

-2.5 

3.77 

6.27 

9.41 

0.75 

3.27 

4.10 

12.00 

0.0 

4.16 

4.16 

8.24 

0.93 

3.53 

4.64 

12.00 

0.0 

4.38 

4.38 

r" . 

8.37 

0.92 

3.64 

4.62 

12.40 

5.35 

2.85 

6.83 

1.20 

4.35 

5.64 

12.40 

2.5 

5.48 

2.98 

6.83 

1.15 

4.16 

5.45 

